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Abstract 

Existence of a new type of oscillating synchronization that oscillates between three different 
types of synchronizations (anticipatory, complete and lag synchronizations) is identified in unidi- 
rectionally coupled nonlinear time-delay systems having two different time-delays, that is feedback 
delay with a periodic delay time modulation and a constant coupling delay. Intermittent anticipa- 
tory, intermittent lag and complete synchronizations are shown to exist in the same system with 
identical delay time modulations in both the delays. The transition from anticipatory to complete 
synchronization and from complete to lag synchronization as a function of coupling delay with 
suitable stability condition is discussed. The intermittent anticipatory and lag synchronizations 
are characterized by the minimum of similarity functions and the intermittent behavior is charac- 
terized by a universal asymptotic — | power law distribution. It is also shown that the delay time 
carved out of the trajectories of the time-delay system with periodic delay time modulation cannot 
be estimated using conventional methods, thereby reducing the possibility of decoding the message 
by phase space reconstruction. 
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Synchronization of chaos is one of the most fundamental phenomena exhib- 
ited by coupled chaotic oscillators. Recent studies on chaotic synchronization 
has also focussed on nonlinear time-delay systems in view of their hyperchaotic 
nature. Further, the concept of delay time modulation has been introduced in 
order to understand dynamical systems with time dependent topology such as 
internet, world wide web, population dynamics, neurology, etc. It has also been 
shown that nonlinear delay systems with time dependent delay can exhibit more 
complex dynamics. Consequently, studies on synchronization of such systems 
with time dependent delay becomes very important in order to understand their 
cooperative dynamics. From this point of view, in this paper we have considered 
simple scalar piecewise linear time-delay systems with unidirectional time-delay 
coupling in order to explore the various types of synchronized behaviours and 
their transitions. The introduction of simple fully rectified sinusoidal modula- 
tion in these systems with constant coupling delay can lead to the existence of 
a new type of oscillating synchronization that oscillates between three diff"erent 
types of synchronizations, namely complete, lag and anticipatory synchroniza- 
tions. With delay time modulation in both the coupling and feedback delay, 
the coupled system displays intermittent lag/anticipatory and complete syn- 
chronizations for suitable ranges of values of the delay times and modulational 
frequencies. It has also been shown that the decrease in the value of modula- 
tional frequency leads to exact and then to lag/anticipatory synchronizations 
from their intermittent nature. The existence of diflFerent synchronizations are 
corroborated by suitable stability condition based on Krasovskii-Lyapunov the- 
ory and their corresponding similarity functions. The intermittent regimes are 
characterized by a universal asymptotic — | power law distribution. 



I. INTRODUCTION 



Chaos synchronization has been receiving a great deal of interest for more than two 
decades in view of its potential applications in various fields of science 
I?]. Since the identification of chaotic synchronization, different kinds of synchronizations 
have been proposed in interacting chaotic systems, which have all been identified both 
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theoretically and experimentally. Complete synchronization refers to a perfect locking of 
chaotic trajectories so as to remain in step with each other in the course of time, X{t) = 
Y{t) 111, 2I, y, I4]. Generalized synchronization is defined as the presence of some functional 



relationship between the states of the drive and the response, Y{t) = F{X{t)) j^, y, 10|. 
Phase synchronization is characterized by entrainment of the phases of the two signals, 
whereas their amplitudes remain uncorrelated ll|, [12] • Lag synchronization implies that 



there is an exact time shift between the evolution of drive and response systems, where the 
response lags the state of the drive, Y{t) = X(t— r) {l3, 14, 15|. Anticipating synchronization 



also appears as a coincidence of shifted-in-time states of t 
anticipates the state of drive, Y(t) = X{t + t 



synchronization has also been identified 



21 



MjlTl, 



18 



le two systems, where the response 
191 . 120| . Recently intermittent lag 



Also it has been shown very recently 



in a three-element network module of semiconductor laser system that dynamical relaying 
can lead to zero-lag synchronization even in the presence of coupling delays 24 1. 

The notion of time dependent delay (TDD) with stochastic or chaotic modulation in time- 
delay systems was introduced by Kye et al 25|] to understand the behaviour of dynamical 
systems with time dependent topology. These authors have reported that in a time-delay 
system with TDD, the reconstructed phase trajectory does not collapse to a simple manifold, 
a property different from that of delayed systems with fixed delay time (which is considered 
to be a serious drawback of the later type of systems). It has been shown very recently 
that a distributed delay enriches the characteristic features of the delayed system over that 
of the fixed delay systems [26]. Based on these considerations, current studies on chaotic 
synchronization in time-delay systems are focused towards time-delay systems with time 



dependent delay 



27 
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301 1 . In this connection it is also of considerable interest to 



study the effect of simple modulations such as periodic modulation |29l . |30| on the nature of 
the chaotic attractor. 

Recently, we have studied chaotic synchronization in a system of two unidirectionally cou- 
pled odd piecewise linear time-delay systems [31] with two different constant delay times: 
one in the coupling term and the other in the individual systems, namely, feedback delay. We 
have shown that there is a transition from anticipatory to lag synchronization through com- 
plete synchronization as a function of a system parameter with suitable stability criterion. 
The present work was motivated by the fact that whether there arises any new phenomenon 
due to the introduction of periodic delay time modulation in the coupled time-delay system 
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we have studied earlier and its effects on the various synchronization scenario. Interestingly, 
we have found that even with simple periodic modulation, the time-delay system cannot be 
collapsed into a simple manifold and that the delay time cannot be extracted using standard 
methods. More interestingly, we have also found that the fully rectified sinusoidal modula- 
tion of delay time introduces a new type of oscillating synchronization that oscillates between 
anticipatory, complete and lag synchronizations for the case of constant coupling delay. This 
is further corroborated by suitable stability condition based on Krasovskii-Lyapunov theory. 
Intermittent anticipatory and lag synchronizations are also found to exist in the present 
system for the case of identical modulation in both the coupling and feedback delay, for a 
range of modulational frequencies. In addition, we also find that there exist regions of ex- 
act anticipatory and lag synchronizations for lower values of modulational frequencies. The 
results have been corroborated by the nature of similarity functions and the intermittent 
behavior by the probability distribution of the laminar phase, satisfying universal — | power 
law behavior of on-off intermittency 32, ISSj . 

The plan of the paper is as follows. In Sec. II, we introduce the scalar piece-wise linear 
time-delay system with delay time modulation and explore the dynamical change in the time 
series of the time-delay system due to delay time modulation. In Sec. Ill we have introduced 
a unidirectional time-delay coupling with delay time modulation between the two scalar 
time-delay systems and we have identified the condition for stability of the synchronized 
state following Krasovskii-Lyapunov theory. A new type of oscillating synchronization that 
oscillates between anticipatory, complete and lag synchronizations and vice versa is shown to 
exist in Sec. IV for the case of constant coupling delay and with modulated feedback delay. 
In Sec. V, we have pointed out the existence of intermittent anticipatory synchronization 
when the strength of the coupling delay is less than that of the feedback delay with identical 
modulations, while in Sec. VI, complete synchronization is realized when the two delays are 
equal. Intermittent lag synchronization is shown to set in when the coupling delay exceeds 
the feedback delay in Sec. VII. In Sec. VIII we very briefly indicate the possibility of more 
complicated oscillating synchronizations in the case of nonidentical modulations. Finally in 
Sec. IX, we summarize our results. 
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II. PIECEWISE LINEAR TIME-DELAY SYSTEM WITH DELAY TIME MODU- 
LATION AND DYNAMICAL CHANGES 



At first, we will introduce the single scalar time-delay system with piecewise linearity 
in the presence of delay time modulation, which has been studied in detail for its chaotic 



dynamics in references 



3, 



35 



36l | with constant time-delay. Then some measures to estimate 



the delay time will be discussed both in the presence and in the absence of delay time 
modulation to show that the imprints of the delay time carved out of the time series of the 
chaotic attractor are completely wiped out by the modulation of delay time. 



A. The scalar delay system 



He 



We consider the following first order delay differential equation introduced by Lu and 
36] and discussed in detail by Thangavel et al. 3^ . 



x{t) = -ax{t) + bf{x{t -t)), 



(1) 



where a and b are parameters, r is the constant time-delay and / is an odd piecewise linear 
function defined as 

0, x< -4/3 

-1.5x - 2, -4/3 <x< -0.8 
f{x) = { X, -0.8 < X < 0.8 (2) 

-1.5X + 2, 0.8<a;<4/3 
0, x> 4/3 

Recently, we have reported [35|] that system ([1]) exhibits hyperchaotic behavior for the pa- 
rameter values a = 1.0, b = 1.2 and r = 25.0 and the hyperchaotic nature was confirmed by 
the existence of multiple positive Lyapunov exponents (see Figures 1 and 2 in ref. (35|). 
Now, we wish to replace the time-delay parameter r as a function of time for our present 



study, instead of the constant time-delay, in the form 



29 



3Q| 



rit) = To + Ta \sm{ujt)\ 



(3) 



where tq is the zero frequency component, is the amplitude and u/tt is the frequency of the 
modulation. Note that in the delay term, we have introduced the fully rectified sinusoidal 



modulational form (absolute of the sine term) so as to keep the delay time positive even 
for values of > tq so as to avoid acausality problem in Eq. ([1]) for negative values of r 
when Ta > Tq. However, for values of Tq sufficiently greater than the rectification in the 
modulation Q is not required. 



B. Estimation of the effect of delay time modulation 



Recently, the concept of time 



was introduced by Kye et al. 



dependent delay with stochastic or chaotic modulation 



25j in the time-delay systems and they have shown in the 
case of Mackey-Glass system that the delay time carved out of time series of the time- 
delay system is undetectable by the conventional measures and hence any reconstruction of 
)hase space of the delayed system is hardly possible. This fact has motivated some authors 

3, y, y, Q ;o look for delay systems with delay time modulation as an ideal candidate 

for secure communication. 

Interestingly we find here that even with a fully rectified sinusoidal delay time modulation 
of the form ([3]), system ([1]) exhibits the properties studied by Kye et al with stochastic or 
chaotic modulation. In order to demonstrate the effect of fully rectified sinusoidal delay 
time modulation of the form ([3]) on the time series of the piecewise linear time-delay system 
which we have considered here, we will calculate (1) filling factor |37l |. (2) length of polygon 
line and (3) average mutual information 25, 39, 40| both in the presence and in the 
absence of delay time modulation and show how periodic modulation removes any imprints 
of the time-delay. 



1. Filling factor 



Now we will compute the filling factor 37|] for the chaotic trajectory x(t) of the time- 
delay system ([1]) by projecting it onto the pseudospace (a;, with P'^^ equally sized 
hypercubes, where the delayed time series Xf = x{t — f) is constructed from x{t) for various 
values of f. The filling factor is the number of hypercubes which are visited by the projected 
trajectory, normalized to the total number of hypercubes, P^^. Figured^ shows the filling 
factor for constant delay tq = 10 when = in Eq. Q, where one can identify the 
existence of an underlying time-delay induced instability |37| which induces local minima 
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in the filling factor at f nro, n = 1,2, 3.... From the later, one can identify the value of 
the time-delay parameter r of the system under consideration. Figure [T]d shows filling 
factor with delay time modulation of the form Q with Tq = 10, = 90 and u = 0.0001, 
where no local minima occurs. Figure [Tt is plotted for tq = 100 and = to show that 
the disappearance of local minima in Fig. [T]d is not due to large delays but only because of 
delay time modulation. From the figures one can realize that the imprints of the delay time 
embedded in the projected trajectory is completely ironed out due to the presence of delay 
time modulation. 



2. Length of polygon line 



Next, to calculate the length of polygon line [38|, the trajectory in {x,Xr,x) space is re- 
stricted to a two dimensional surface. The restriction in dimension is effected by intersecting 
the projected trajectory with a surface k{x,Xr,x) = 0. Consequently the number of times 
the trajectory traverses the surface and the corresponding intersection points can be calcu- 
lated. One then orders the points with respect to the values of , and a simple measure for 
the alignment of the points is the length L of polygon line connecting all the ordered points. 
Figure [2^ shows length of polygon line L with constant delay Tq = 10, where the local min- 
ima correspond to the delay time of the system we have considered. Figure [2)3 shows length 
of polygon line L with delay time modulation where there is no remnance of information 
about delay time from the trajectory, whereas Fig. [2t is plotted for tq = 100, = 0, to show 
that the imprints of delay time carved out in the trajectory vanishes in Fig. [2]d only due to 
the delay time modulation and not because of large delay. 

3. Average mutual information 



As a final example, we will calculate average mutual information defined by (see for 



example. 



25 



39 



40(1 and references therein) 



H^) = Yl Pi^i'"'), xin + f)) X log2 



P{x{n), x{n + f)) 
P{x{n))P{x{n + f)) 



(4) 



where P{x{n),x{n + f)) is the joint probability density for measurements in the chaotic 
time series X = (a;(l), a;(2), ...,x{m)) and in the constructed delay time series = (x(l + 
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f),x{2 + f), ...,x{m + f)) by varying f, resulting in values x{n) and x{n + f). P{x{n)) and 
P{x{n + f)) are the individual probability densities for the measurements of X and Xf. 
Figure [3] shows the average mutual information for the cases of constant delay time with 
tq = 10 (Fig. [3^1') and with delay time modulation (Fig. [nb)- Figure|3t is plotted for tq = 100 
to show that the absence of local peaks in Fig. [3)3 is due to delay time modulation and not 
because of large delay. For fixed delay time the average mutual information shows local 
peaks at the time-delay f = tq (or multiples of it f = nro) of the system, whereas for the 
case of delay time modulation the average mutual information has no such peaks to identify 
the delay time of the delayed system. 

One can also obtain similar results with other measures such as autocorrelation func- 
tion, onestep prediction error and average fitting error 



33, 



38 



41 



42| . However, we are 



not presenting these results here for convenience. In order to perform the phase space re- 
construction, the first step is to find out the delay time for the projected trajectories. By 
introducing the delay time modulation the imprints of delay time in the projected trajec- 
tory is completely removed as seen above for the present system, inhibiting any possibility 
of phase space reconstruction. This is essentially consequent of the fact that when the delay 
time is modulated by the fully rectified sine term, the delay time effectively gets increased 
in which case the number of positive Lyapunov exponent also increases (as noted in Fig. 2 



in Ref. 



35|. Consequently study of chaos synchronization in a system of such coupled delay 



time modulated oscillators will be of considerable interest. 



III. COUPLED SYSTEM AND THE STABILITY CONDITION IN THE PRES- 
ENCE OF DELAY TIME MODULATION 

Now let us consider the following unidirectionally coupled drive Xi{t) and response X2{t) 
systems with two different modulated time-delays Ti{t) and T2{t) as feedback and coupling 
time-delays, respectively (hereafter we write Ti(t) and T2(t) simply as ri and T2 respectively), 

xi{t) = -axi{t) + bif{xi{t - Ti)), (5a) 
X2{t) = -ax2{t) + h2f{x2{t - n)) + hf{xi{t - rs)), (5b) 
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where bi, 62 and 63 are constants, a > 0, and f{x) is of the same form as in Eq. ([2]) with 

n = Tio + \sm{ujit) \ , (6a) 
T2 = r2o + |sin(u;2t)| , (6b) 

where tiq and T20 are the zero frequency components of feedback delay and couphng delay, 
Tia and T2a are the amplitudes of the time dependent components of ri and T2, respectively, 
and cui/vr and 002/ tt are the corresponding frequencies of their modulations. 

Now we can deduce the stability condition for synchronization of the two time-delay 
systems, Eqs. (EK) and (Eb), in the presence of the delay coupling h^f{xi{t — T2)) with time 
delay modulation in both the feedback delay and coupling delay. The time evolution of the 
difference system with the state variable A = Xir^^r^ — ^2, where Xi^-a-n = Xi{t — — ti)) , 
can be written for small values of A by using the evolution Eqs. ([5]) as 

A = -aA + (62 + h- bi)f{xi{t - r2)) + b2f'{xi{t - r2))A,,, A. = A{t - r). (7) 

Then A = corresponds to anticipatory synchronization when T2 < ti, identical or complete 
synchronization for T2 = Ti and lag synchronization when T2 > Ti. In order to study the 
stability of the synchronization manifold as in the case of constant time delay case [sil, we 
choose the parametric condition, 

61 = 62 + &3, (8) 
so that the evolution equation for the difference system A becomes 

A = -aA + 62/'(xi(t-r2))A.,. (9) 
The synchronization manifold is locally attracting if the origin of this equation is stable. 



43 



4J], we define a positive definite 



Following Krasovskii-Lyapunov functional approach 
Lyapunov functional of the form 

V{t) = ]-A^ + f A\t + 9)d9, (10) 

where /i is an arbitrary positive parameter, yU > 0. Note that V{t) approaches zero as A ^ 0. 

To estimate a sufficient condition for the stability of the solution A = 0, we require the 
derivative of the functional V(t) along the trajectory of Eq. ([9]), 

^ = -aA' + b2f'Mt - r2))AA., + [Alj[ + A' - A^J , (11) 
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to be negative. Note that in the case of constant modulation t[ = ^ vanishes. The above 
equation can be rewritten as 



(12) 



wherer =[((a-/i)//i)-(62/'(a;i(t-r2))//i)X + XV(l-r{)], X = A,,/A. In order to show 
that ^ < for all A and A^- and so for all X, it is sufficient to show that T^in > 0. One 
can easily check that the absolute minimum of F occurs at X = b2f'{xi{t — r2))/2/i(l — t[) 
with Tmin =[4/i(a — yu)(l — r[) — blf''^{xi{t — r2))]/4yU,^(l — t[). Consequently, we have the 
condition for stability as 



a > 



4/i(l 



(13) 



Again $(/i) as a function of fi for a given f'{x) has an absolute minimum at /i 



b2f'{xi{t-T2)) 



with <l>r 



b2f'ixi{t-T2)) 



Since $ > <l>r 



b2f'{xi{t-T2)) 



from the 



inequahty (fT3l) . it turns out that the sufficient condition for asymptotic stabihty is 



a > 



(14) 



along with the condition ([8]) on the parameters bi, 62 and 63. 

Now from the form of the piecewise linear function f{x) given by Eq. (2), we have, 



\f{x^{t-T,))\ 



1.5, 0.8 < < I 



1.0, 



\xi\ < 0.8 



(15) 



Consequently the stability condition (IT^ becomes a > 1.5 



the parametric restriction bi 



Thus one can take a > 



b2 



b2 



> 



b2 



along with 



as a less stringent condition for (HM to be valid, while 



a > 1.5 



(16) 



can be considered as the most general condition specified by (1141) for asymptotic stability 
of the synchronized state A = 0. The condition ( |T6i) indeed corresponds to the stability 
condition for exact anticipatory/lag as well as exact complete synchronizations for a given 
value of the coupling delay T2 in a global sense. It may be noted that the stability condition 
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(fT6|) is valid irrespective of the nature of the couphng delay, that is whether it is constant or 
modulated. However, when the feed back delay ri is constant the condition f|T6l) reduces to 
a > 1.5 1 62 1 as discussed in ref. 311- In the following, we will consider both the possibilities of 
constant {t2 = T20) and periodically modulated (r2 = T2o+T2a |sin(ci;2^)|) coupling delays with 
a periodically modulated feedback delay (ri = tiq + Tia |sin(c<Jit)|). We demonstrate through 
detailed numerical analysis that there exists oscillating synchronization that oscillates be- 
tween anticipatory, complete and lag synchronizations for the case of constant coupling delay 
T2 = '^20- Intermittent anticipatory/lag and complete synchronizations are shown to exist 
for the case of coupling delay with delay time modulation T2 = T20 + T2a |sin(ci;2t)|, when 
T2a = Tia and oji = ijJ2- For r2a 7^ T\a and uj\ 7^ u;2, more complicated oscillating type 
synchronizations occur. 



IV. OSCILLATING SYNCHRONIZATION 



At first we consider the constant coupling delay, T2 = T20, and show that there exists 
oscillating synchronization that oscillates between anticipatory, complete and lag synchro- 
nizations as a function of time for suitable range of parameters. 

Now we will choose the delay time modulation in the form ([6^) for the feedback delay ri(= 
Tio + ^ia |sin(u;it)|) with tiq = 10, Tia = 90 and Ui = 10~'^. We have fixed the value of T2a = 
in (Eb), so that the coupling delay becomes constant T2 = T20 = 45 with the parameters 
a = 1, 61 = 1.2 in Eq. ([5]) and the values of 62 and 63 are chosen according to the parametric 
restriction ([8]) depending upon the stability condition to be satisfied. For the chosen values 
of Tio and Tia, one can find that Ti oscillates between (ri(t) = Tiq + Tia |sin(ct;it)| = 10 + 
90 |sin(u;it)|) the values 10 and 100. With the chosen value of constant coupling delay 
T2 = 45 and time dependent feedback delay ti, as time evolves one finds that the feedback 
delay Ti(t) is lesser than the value of constant coupling delay T2 initially for some time (in 
which case r(t) = T2 — Ti(t) > 0, so that there exists lag synchronization xi{t — T{t)) = X2{t) 
with varying lag time r(t) = T2 — Ti{t)). As time evolves, Ti{t) increases eventually and it 
approaches ti = 45 at a certain later time (T = n/uJi), where r(t) = T2 — Ti{t) = 0, so 
that Xi{t) = X2{t) and a complete synchronization occurs at a specific value of time. As 
Ti{t) increases further above the value of T2 = 45, the delay time r(t) becomes negative, 
r(t) = T2 — Ti{t) < with xi{t — r{t)) = X2{t)) and there exists anticipatory synchronization 
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with varying anticipating time r(t) = T2 — Ti{t). This anticipatory synchronization continues 
till the value of time dependent feed back delay ri (t) decreases to appraoch the value of the 
constant coupling delay T2 = 10 after reaching its maximum value of 100. Therefore as time 
evolves there is oscillation between lag, complete and anticipatory synchronizations with 
time dependent anticipating and lag times. 

Figure UK shows the evolution of the drive Xi{t) and the response X2{t) at the transi- 
tion between lag to anticipatory synchronization via complete synchronization for the value 
of 62 = 0.1, where the general stability condition (fT6|) is satisfied, whereas Fig. |Dd shows 
the evolution of the drive Xi{t) and the response X2{t) at the next transition between an- 
ticipatory to lag via complete synchronization. In Figs. and \5jp, the difference signals 
Xi{t — t) — X2(t) , T > and Xi(t — r) — X2(t) , r < are plotted respectively for the value of 
parameters satisfying the general stability condition corresponding to the Fig. |U confirm- 
ing the transition between lag to anticipatory synchronization. Thus as a consequence of 
delay time modulation there exists a new type of oscillating synchronization that oscillates 
between anticipatory, complete and lag synchronizations with varying anticipating and lag 
times. 

V. INTERMITTENT ANTICIPATORY SYNCHRONIZATION 

Now we consider the coupled time-delay system ([5]) with delay time modulation of the 
form in both the feedback and coupling delays for further studies. We have fixed the 
values of the parameters as a = l,6i = 1.2, = T2a = 90, Ui = uj2 = 10~^ (identical 
modulations) and the values of 62 and 63 are chosen according to the parametric restriction 
bi = 62 + ^3 depending upon the stability condition to be satisfied. For ri, the zero frequency 
component of amplitude is fixed as tiq = 10 and for T2, it is fixed as T20 = 5, so that a constant 
difference is maintained between the feedback and the coupling time delays throughout the 
time evolution. With the coupling delay T2{= 5+90 |sin(10~^t)|) being less than the feedback 
delay ti(= 10 + 90 |sin(10~^t)|), that is r2(t) < Ti(t), the value of the anticipating time r = 
T2 — Ti turns out to be negative such that the relation between drive xi{t) and the response 
X2{t) now becomes xi{t — r) = X2{t),T < 0, demonstrating anticipatory synchronization, 
provided the stability condition flTB]) is satisfied with the parametric restriction specified by 
Eq. (ED. 
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Now let us choose the parameter 62 as the control parameter, whose value determines the 



1. For 62 = 
with \f\ 



0.7, 1.5 



by Eq. 








> a > 









the less stringent condition is satisfied 



-( ^ 1 for the chosen values of u and r^. One can observe intermittent 
anticipatory synchronization as shown in Fig. [6l exhibiting typical features of on-off 



intermittency 32, l33| with the off state near the laminar phase and the on state 
showing a random burst. For this value of 62 the amplitude of the laminar phase 
corresponding to the synchronized state is approximately zero (of the order 10~^). 



2. For bo 



0.1, a > 1.5 



62 



> 



b2 



the general stability condition (1161) is sat- 
analysis reveals that here the intermittent 



isfied and correspondingly the numerica 
anticipatory synchronization is such that the amplitude of the laminar phases corre- 
sponding to the synchronized state is exactly zero (in the sense that the difference 
A = xi(t — r) — X2(t), T < is of the order 10~^^ in the laminar phases) as shown in 
Fig. El 

To analyze the statistical features associated with the intermittent nature in Fig. [7] for 
the value of 62 = 0.1, we have calculated the distribution of laminar phases A{t) with 
the amplitude less than a threshold value A < lO^^'' and we have observed a universal 



asymptotic 



power law distribution as shown in Fig. [HI which is quite typical for on-off 



32 



33| . One can also find a similar power law distribution for the value of 



mtermittency 

62 = 0.7 discussed above but now with a larger threshold value (A < 10~^) of the laminar 
region. 



Now using the notion of similarity function introduced by Rosenblum et al. HI] to 
characterize lag synchronization, one can also characterize anticipatory synchronization. 
Similarity function for anticipatory synchronization is defined as the time-averaged difference 
between the variables Xi and X2 (with mean values being subtracted) taken with the time 
shift |t|, 

([xi(t+|r|)-X2(t)n 



Slir) 



(17) 



where, (x) means time average over the variable x. If the minimum value of (r) reaches 
zero, that is S'a(r) = 0, then there exists a time shift |r| between the two signals Xi{t) and 



13 



X2{t) such that xi{t+ |r|) = X2{t), demonstrating the existence of anticipatory synchroniza- 
tion between the drive xi and the response X2 signals. Figure [9] shows the similarity function 
S'a(r) as a function of the difference between the feedback and the coupling delays, r = T2—T1 
for three different values of 62? the parameter whose value determines the stability condition 



b2 



> 



b2 



> a) 



In Fig. [HI the Curve 3 is plotted for the value of 62 = 1-1, (1-5 
where both the less stringent condition and the most general condition are violated. From 
the curve 3 one can find that the minimum value of 5'a(r) is greater than zero for all values 
of r, resulting in the lack of exact time shift (anticipating time) between the drive and 
the response signals. On the other hand the curve 2 corresponds to the value of 62 = 0.7 
such that the less stringent condition is satisfied while the general stability criterion fll6p 
is violated as seen above. Curve 2 shows that the minimum of similarity function ^^(r) is 
approximately zero (of the order 10"'^) for r < 0, as may be seen in the inset of Fig. [9l 
indicating the existence of intermittent anticipatory synchronization with the amplitude of 
the laminar phases of the difference signal A = xi{t — t) — X2{t) , t < 0, being approximately 
zero (< 10^^). On the other hand, the curve 1 is plotted for 62 = 0.1, satisfying the general 
stability criterion Eq. ( |T6l) . which shows that the minimum of similarity function is much 
closer to zero (of order 10~^), r < 0, indicating that there exists an intermittent anticipatory 
synchronization with the amplitude of the laminar phase of the difference signal becoming 
exactly zero with the anticipating time equal to the difference between the two time delays 

T = T2- Ti. 

Next, by reducing the value of the modulational frequencies u = ui = 002 further, we find 
that the lengths of the laminar phases increase gradually with a corresponding decrease in the 
number of turbulent phases. Finally at an appropriate value of the modulational frequency 
all the turbulent phases disappear and there exists only exact anticipatory synchronization 
without any intermittent bursts provided the exact stability condition is satisfied. Corre- 
spondingly the similarity function S'„,(r) becomes zero exactly (which is of the order 10~^^) 
for r < in this case, as shown in [sij. 

VI. COMPLETE SYNCHRONIZATION 

Complete synchronization follows the anticipatory synchronization when the value of the 
coupling time-delay T2 equals the feedback time-delay ri, that is T2 = ri, where the antici- 
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pating time becomes t = T2 — ti = 0. Here also, the same stability criterion Eq. (fT6!) holds 
good with the same parametric restriction specified by ([8]). In this case of complete synchro- 
nization {t2 = Ti), the delay time modulation does not induce any intermittent nature in 
the dynamical behavior of the coupled systems for any value of the modulational frequency 
(cui = UJ2) as inferred from Eq. (6). Figure [TUk shows as an illustration the plot of Xi(t) vs 
X2{t) for the values of 62 = 0.7 and ui = UJ2 = 10~^, such that the less stringent condition is 
satisfied and the general stability criterion f|T6l) is violated. The plot shows small deviations 
from the localized diagonal line implying an approximate synchronization, whereas Fig. [TOb 
shows an entirely localized sharp diagonal line for the value of 62 = 0.1, where the general 
stability condition (ITBi) is satisfied, indicating the complete synchronization. 

VII. INTERMITTENT LAG SYNCHRONIZATION 

When the value of the coupling delay T2 is increased above the value of the feedback delay 
ti(t2 > Ti), then the value of the retarded time t = r2 — ti turns out to be positive such 
that the relation between the drive Xi(t) and the response X2(t) now becomes Xi(t — r) = 
X2{t),T > 0, depicting the existence of lag synchronization, provided the general stability 
condition (fT6!) is satisfied along with the parametric condition ([8]). 

We have fixed the same values for all the parameters as in the case of intermittent 
anticipatory synchronization except for the zero frequency component T20 of coupling delay 
T2 which is fixed at T20 = 15. Figure [TT] shows the intermittent lag synchronization for the 
value of 62 = 0.7, in which case only the less stringent stability condition is satisfied, where 
the laminar phase has an amplitude which is nearly zero (of the order 10^^). Figure [T2] 
shows intermittent lag synchronization for the value of 62 = 0.1, where the amplitude of the 
laminar phase vanishes exactly. In the later case the most general stability criterion (ITB]) 
is satisfied. The statistical behavior associated with the intermittent nature in this case 
of intermittent lag synchronization is also characterized by the probability distribution of 
laminar phases having amplitude less than a threshold value A < 10^^° corresponding to a 
universal asymptotic — | power law distribution as shown in the Fig. [131 

The figure shows the probability distribution A(t) of intermittent lag synchronization 
for the value of 62 = 0.1. One can also verify that the intermittent lag synchronization 
for the value of 62 = 0.7 has also similar power law distribution for larger threshold value 
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(A < 10~^)of amplitude of the laminar phases. 

The existence of intermittent lag synchronization is also characterized by a similarity 
function Si{t) defined as 

N _ ([xi(t- |r|) -X2(t)]') 

Figure UM shows the similarity function Si{t) for intermittent lag synchronization as a func- 
tion of the retarded time t = T2 — ti. Curve 3 is plotted for the value of 62 = 1-1 (which is 



greater than both ayl — t[ and ay^l — t[/1.5), where the minimum of similarity function 
Si{t) occurs at a finite value of Si{t) > and hence there is a lack of lag synchronization 
between the drive and the response signals indicating asynchronization. Curve 2 corresponds 



to the value of 62 = 0.7, (which is less than a^yl — t[ but greater than a^/l — r{/1.5), where 
the minimum of similarity function Si{t) is approximately zero (of the order of 10~^, as may 
be seen in the inset of Fig. [T^ indicating the existence of intermittent lag synchronization 
with the amplitude of the laminar phase being approximately zero. However, for the value of 
62 = 0.1, for which the general condition f|T6|) is obeyed, the minimum of similarity function 
for Curve 1 becomes much closer to zero (of the order 10^^) which corresponds to intermit- 
tent lag synchronization with exact time shift between the two signals during the laminar 
phase. 

Next, as in the case of intermittent anticipatory synchronization, by reducing the value 
of modulational frequency one can find that the lengths of the laminar phases increase with 
vanishing turbulent phases and finally at an appropriate value of the modulational frequency 
there exists exact lag synchronization without any intermittent bursts provided the exact 
stability condition is satisfied. Correspondingly the similarity function Si{t) becomes zero 
exactly (which is of the order 10~^^) for r > in this case. 



VIII. COMPLEX OSCILLATING SYNCHRONIZATION 

Finally, when 7^ or /and Ui 7^ uj2 the frequencies as well as amplitudes of the 
modulated feedback delay ri(t)(= tiq + Tia\sm{uit)\) and the modulated coupling delay 
''"2(^)(= T20 + '^2a |sin(ci;2^)|) differ from each other resulting in a more complicated variation 
of the anticipating/lag time T{t) = T2(t) — Ti{t). This in turn results in the existence of 
more complex oscillating synchronization than the one presented in Sec. IV. It is clear that 
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one can also introduce other kinds of modulations instead of periodic modulation to obtain 
varying forms of oscillating synchronizations. 

IX. SUMMARY AND CONCLUSIONS 

In this paper, we have shown that there exists a new type of oscillating synchronization 
that oscillates between anticipatory, complete and lag synchronization and vice versa for 
the case of constant coupling delay with varying anticipating and lag times. We have also 
shown that there exists regions of intermittent anticipatory synchronization, intermittent 
lag synchronization and complete synchronization in the parameter space of uj and T2 with 
appropriate stability condition for the synchronized state in a system of two piecewise linear 
time-delay systems with delay time modulation in both the feedback and coupling delay. For 
a fixed value of tu, we have shown that there is a transition from intermittent anticipatory to 
intermittent lag synchronization through complete synchronization with the coupling delay 
T2 as the only control parameter, while all the other parameters are kept fixed. We have 
also found that on further reducing the value of tu, one can observe transition towards ex- 
act anticipatory/lag synchronizations from their intermittent behaviour. The signature of 
the intermittent behavior in both the intermittent anticipatory and intermittent lag syn- 
chronizations are characterized by probability distribution of laminar phases satisfying a 
universal asymptotic — | power law distribution. The existence of intermittent anticipatory 
and intermittent lag synchronizations are characterized by their corresponding similarity 
functions. 

Further, we have observed that in the region where the stringent stability condition (fT6|) 
is satisfied, the minimum of the similarity functions S'a(r) and Si{t) approaches very closely 
zero for all values of T2 < ti and T2 > ri, respectively. The range of zero values corresponding 
to the minimum of similarity functions S'a(r) and Si{t) shows the existence of anticipatory 
and lag synchronizations for the values of coupling delay T2 below and above the feedback 
delay ri, respectively. We have also shown that the estimation of the delay time carved out 
of the time series of the delayed system even with delay time modulations of fully rectified 
sinusoidal type is very difficult by conventional methods in the present system for suitable 
choice of the parameters (in contrast to the chaotic or stochastic delay time modulation as 
studied by Kye et al [25] in Mackey-Glass delay system) and so the messages encoded in such 
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systems can be expected to be less amenable for extraction by phase space reconstruction. 
We have also confirmed numerically that the phenomena reported in this paper occur in 
other time-delay systems such as Mackey-Glass and Ikeda systems also to corroborate the 
generic nature of the results. Also the model system discussed in the present manuscript 
is amenable for experimental realization in terms of suitable nonlinear electronic circuits in 
view of its piecewise linear nature and we are pursuing the experimental verification of it. 
It is hoped that the study of such modulated systems will lead to a better understanding 
of the dynamics of systems with time dependent topologies such as neural networks, world 
wide web, etc. 
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Figure captions 



1. Filling factor as a function of delay time f (of delayed time series Xr). (a) with 
constant delay tq = 10 when = 0, (b) with delay time modulation of the form 
with To = 10, Ta = 90 and u = 10"^ and (c) with large constant delay tq = 100(ra = 0). 

2. Length of polygon line as a function of delay time f (of delayed time series Xf ). (a) with 
constant delay tq = 10 (r^ = 0), (b) with delay time modulation of the form with 
the same parameters as in Fig. [T]and (c) with large constant delay tq = 100 (r^ = 0). 

3. Average mutual information as a function of delay time f (of delayed time series 
Xr). (a) with constant delay tq = 10, = 0, (b) with delay time modulation of the 
form ([3]) with the same parameters as in Fig. [T] and (c) with large constant delay 

To = 100, Ta = 0. 

4. Oscillating synchronization for the constant coupling delay T2 = 45 with time depen- 
dent feedback delay of the form (jH^) with tiq = 10, Tia = 90 and u = 10^^. (a) 
Oscillating from lag to anticipatory synchronization via complete synchronization in 
the region t G (3970, 4020) and (b) Oscillating from anticipatory to lag synchronization 
at the next transition in the region t G (27400, 27450). 

5. (a) Difference between Xi(t — r),r > and X2{t), showing lag synchronization for 
certain time and (b) difference between a;i(t — r), r < and X2{t), showing anticipatory 
synchronization for the following period of time for 62 = 0.1 satisfying the general 
stability condition (flGll . Note that complete synchronization occurs in the transition 
regime. 

6. The time series Xi{t — r) — X2{t),T < 0, for 62 = 0.7 and 63 = 0.5 (so that the less 
stringent condition a > I62/ \/l ~ ^il is satisfied while (fT6ll is violated) corresponding 
to intermittent anticipatory synchronization with the amplitude of the laminar phase 
approximately zero. 

7. The time series Xi{t — r) — X2{t),T < 0, for 62 = 0.1 and 63 = 1.1. Here the general 
stability criterion (fT6l) is satisfied corresponding to intermittent anticipatory synchro- 
nization with the amplitude of the laminar phase exactly zero. 
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8. The statistical distribution of laminar phase satisfying — | power law scaling for 62 = 
0.1 and 63 = 1.1, where the general stability criterion f[T6l) is satisfied. 

9. Similarity function for intermittent anticipatory synchronization S'a(r) for different 
values of 62? the other system parameters are a = 1.0, 61 = 1.2 and u = 10^^. (Curve 1: 
62 = 0.1,63 = 1.1, Curve 2: 62 = 0.7,63 = 0.5 and Curve 3: 62 = 1.1,63 = 0.1). 

10. Complete synchronization between Xi{t) vs X2{t) when T20 = tiq. (a) Approximate 
complete synchronization for 62 = 0.7 and (b) Exact complete synchronization for 
62 = 0.1. 

11. The time series Xi{t — r) — 0:2(^)7''" > 0, for 62 = 0.7 and 63 = 0.5 (so that the less 
stringent condition a > [62/ ^Jl — t[\ is satisfied while f|T6|) is violated) corresponding 
to intermittent lag synchronization with the amplitude of the laminar phase approxi- 
mately zero. 

12. The time series xi{t — r) — X2{t),T > 0, for 62 = 0.1 and 63 = 1.1. Here the general 
stability criterion f|T6l) is satisfied corresponding to intermittent lag synchronization 
with the amplitude of the laminar phase exactly zero. 

13. The statistical distribution of laminar phase satisfying — | power law scaling for 62 = 
0.1 and 63 = 1.1, where the general stability criterion (fT6l) is satisfied. 

14. Similarity function for intermittent lag synchronization Si{t) for different values of 
62, the other system parameters are a = 1.0, 61 = 1.2 and u = 10^^. (Curve 1: 
62 = 0.1, 63 = 1.1, Curve 2: 62 = 0.7, 63 = 0.5 and Curve 3: 62 = 1.1, 63 = 1.0). 
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